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Abstract 

\Q . 

. In this paper we explore the possibility of suitably combining the 

CNj ' nodes SI of the existing geodetic LAGEOS, LAGEOS II and Ajisai 

^^O . laser-ranged satellites and of the radar altimeter Jason-1 satellite in 

order to increase the accuracy in testing the general relativistic gravit- 
omagnetic Lense-Thirring secular effect in the gravitational field of the 
' Earth. The proposal of introducing Ajisai and Jason-1 in such a combi- 

nation comes from the expected benefits which could be obtained in re- 
' ducing the aliasing secular impact of the classical part of the terrestrial 

Qi. gravitational field. According to the recently released EIGEN-CG01C 

^ ' combined GRACE+CHAMP+terrestrial gravimetry/altimetry Earth 

bi)' gravity model, the impact of the static part of the mismodelled even 

^ ' zonal harmonics of geopotential, which represent the major source of 

•i— i . systematic error, amounts to 1.6%, at 1— a level. It is better than the 

' error which could be obtained with a two-node LAGEOS-LAGEOS 

II only combination (6% at 1— <r). Moreover, the proposed combi- 
nation would be insensitive also to the secular variations of the low- 
degree even zonal harmonics, contrary to the LAGEOS-LAGEOS II 
only combination. Such variations could be a serious limiting factor 
over observational time spans many years long. The price to be paid 
for this improvements of the systematic error of gravitational origin is 
represented by the non-conservative forces introduced along with the 
new orbital elements. However, they would induce periodic pertur- 
bations, contrary to the gravitational noise. A major concern would 
be the assessment of the impact of the non-conservative accelerations 
on the Jason-1 node. According to the present-day force models, the 
mismodelling in the non-conservative forces would, at worst, induce an 
aliasing periodic signal with an amplitude of 4% of the Lense-Thirring 
effect over a time span of 2 years. However, an observational time 
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span of just some years could safely be adopted in order to fit and re- 
move the residual long-period non-gravitational signals affecting Ja- 
son's node, which, in the case of the direct solar radiation pressure, 
have a main periodicity of approximately 120 days. Of course, the 
possibility of getting time series of the Jason's node some years long 
should be demonstrated in reality. 



1 Introduction 

1.1 The performed attempts to measure the Lense— Thirring 
effect with the LAGEOS satellites 

The general relativistic gravitomagnetic force pQ, related to the component 
of the gravitational field induced by the rotation of a central body of mass 
M and proper angular momentum J, is still awaiting for a direct, unques- 
tionable measurement. Up to now there exist some indirect evidences of 
its existence as predicted by the General Theory of Relativity (GTR in the 
following) in an astrophysical, strong-field context [2 ; and, in the weak-field 
and slow-motion approximation valid throughout the Solar System, in the 
fitting of the ranging data to the orbit of Moon with the Lunar Laser Rang- 
ing (LLR) technique [3]. 

The direct measurement of the gravitomagnetic Schiff precession of the 
spins of four spaceborne gyroscopes [I] in the gravitational field of Earth is 
the goal of the Stanford GP-B mission ,5_ which has been launched in April 
2004. The claimed obtainable accuracy is of the order of 1% or better. 

Another consequence of the gravitomagnetic field is the Lense-Thirring 
effect on the geodesic path of a test particle freely orbiting a central rotat- 
ing body [Hj. It consists of tiny secular precessions of the longitude of the 
ascending node O and the argument of pericentre uj of the orbit of the test 
particle 

■ 2GJ . 6GJcosi 

"lt = r, ^lt = r, (1) 

c 2 a 3 (l — e 2 )a c 2 a 3 (l — e 2 )z 

where a, e and i are the semimajor axis, the eccentricity and the inclination, 
respectively, of the orbit, c is the speed of light and G is the Newtonian grav- 
itational constant. See Figure fTTTl for the orbital geometry. The LAGEOS 
III/LARES mission [7||H] was specifically designed in order to measure such 
effect, but, up to now, in spite of its scientific validity and relatively low 
cost, it has not yet been approved by any space agency or scientific insti- 
tution. Recently, a drag- free version of this project 0, in the context of 
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Figure 1: Orbital geometry for a motion around a central mass. Here L 
denotes the orbital angular momentum of the particle of mass m, J is the 
proper angular momentum of the central mass M, II denotes the pericentre 
position, / is the true anomaly of m, which is counted from II, f2, to and 
i are the longitude of the ascending node, the argument of pericentre and 
the inclination of the orbit with respect to the inertial frame {x, y, z} and 
the azimuthal angle <f> is the right ascension counted from the x axis. When 
orbits with small inclinations are considered, the longitude of pericentre 
w = ft + ijj is used. 
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the relativistic OPTIS mission ^U], is currently under examination by the 
German Space Agency (DLR). 

2 The current LAGEOS-LAGEOS II Lense-Thirring 
experiment 

The Ciufolini's team analyzed the laser-ranged data to the existing geodetic 
LAGEOS and LAGEOS II satellites over time spans of some years |11M12| by 
using the following J2 — J a— free 1 linear combination of the orbital residuals 
of the nodes of LAGEOS and LAGEOS II and the perigee of LAGEOS II 
PI 

Sn hAGEOS + Cl ^ LAGEOS 11 + c 2 ^ LAGEOS 11 -^60.2, (2) 

where c\ = 0.304, C2 = —0.350 and /xlt is the solved-for least squares pa- 
rameter which is in Newtonian mechanics and 1 in the GTR 2 . The resid- 
uals of the nodal and perigee rates 5£l and dui are determined by comparing 
these orbital elements at the end of two orbit computation arcs covering the 
same time span. These orbit computations consist of a precise orbit, which 
has been adjusted to fit the available tracking data, and a predicted orbit, 
starting from an initial state vector which it shares with the precise orbit, 
but for which the later states are based on integration of the classical force 
models, without adjustment to the tracking data. The residuals of the rate 
of change of the orbital elements are then calculated by dividing the differ- 
ence between these two elements by the time span of the arc. The signature 
of the gravitomagnetic force, assumed as an unmodelled effect, is a linear 
trend in the combination of the accumulated residuals, with a slope of 60.2 
milliarcseconds per year (mas yr -1 in the following). 

The latest, 2002, measurement of the Lense-Thirring effect with Pj). 
obtained by processing the LAGEOS and LAGEOS II data over a time span 
of almost 8 years with the orbital processor GEODYN II of the Goddard 
Space Flight Center, yields |T2] 

m ~l±0.02±^ tematic , (3) 

1 J2 and J4 are the first two even zonal harmonic coefficients of the multipolar expansion 
of the terrestrial gravitational field, ft turns out that they induce secular precessions on 
Q and u. Their mismodelled parts, according to the currently available Earth gravity 
models, are of the same order, or even larger, than the gravitomagnetic secular precessions 
of interest. 

2 It can be expressed in terms of the PPN parameter 7 as /ilt = ^p-- For the PPN 
formalism see |14l HI [To|. 
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where <5/i L ^ tematlc accounts for all the possible systematic errors due to the 
mismodelling in the various competing classical forces of gravitational and 
non-gravitational origin affecting the motion of the LAGEOS satellites. In 
[T2] <5//LT tematlC is estimated to be of the order of 



^systematic = ^ _ 3Q%> (4) 



2.1 Some possible criticisms to the performed attempts to 
measure the Lense— Thirring effect with the LAGEOS 
satellites 

However, it must be pointed out that several remarks have been made by 
other scientists about the estimate (JIJ (EE E|- First of all, it is based 
on assuming a 13% systematic error due to the mismodelling in the even 
zonal harmonics of geopotential. This result has been obtained from the full 
covariance matrix of the EGM96 Earth gravity model [IB, up to degree i = 
20 3 . In the EGM96 solution the recovered even zonal harmonics are strongly 
reciprocally correlated; it seems, e.g., that this value for the systematic error 
due to geopotential is due to a lucky correlation between Jq and Js which 
are not cancelled by (J2J). The point is that, according to ^3 Eh nothing 
would assure that the covariance matrix of EGM96, which is based on a 
multi-year average that spans the 1970, 1980 and early 1990 decades, would 
reflect the true correlations between the even zonal harmonics during the 
particular time intervals of a few years adopted in the analyses by Ciufolini 
and coworkers. Then, a more conservative, although pessimistic, approach 
would be to consider the sum of the absolute values of the errors due to 
the single even zonal as representative of the systematic error induced by 
our uncertainty in the terrestrial gravitational field according to EGM96 
^21120]. In this case we would get a pessimistic upper bound of 83%. 

Also the use of the perigee of LAGEOS II in © sheds shadows on the 
estimate of @. Indeed, it turns out that the perigee of the LAGEOS-like 
satellites is very sensitive to a whole host of non-gravitational perturba- 
tions of radiative |2I] and thermal origin j^j , whose impact on the proposed 
measurement of the Lense-Thirring effect would be very difficult to reliably 
assess |161 117j . It should also be pointed out that at the time of the analy- 
ses by Ciufolini and coworkers, non-gravitational perturbations of thermal 
origin, such as the solar Yarkovsky-Schach effect, were not present at all in 

3 In fact, it turns out that the result does not change if the even zonal harmonics of 
degree higher than £ — 12 are neglected in the calculation. 
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the force models of the orbital determination software GEODYN II, which 
was used in their analyses. 

3 The new perspectives opened up by the CHAMP 
and GRACE models 

The recent and forthcoming improvement in our knowledge of the Earth 
gravity field thanks to CHAMP |23] and especially GRACE [T£|, opens 
new possibilities to obtain a more accurate and reliable measurement of 
the Lense-Thirring effect with the LAGEOS satellites as well as the other 
currently existing laser-tracked satellites. In the pre-CHAMP and GRACE 
era these other satellites would have been unsuitable j^l] because of their 
sensitivity to the higher degree even zonal harmonics of geopotential, due 
to their smaller altitude with respect to LAGEOS and LAGEOS II. In 
a new observable has been explicitly put forth, which consists of a J2— free 
combination of the nodes of LAGEOS and LAGEOS II only. It is 4 

^LAGEOS + ^LAGEOS II „ (5) 

with qi = 0.546. According to the recently released EIGEN-CG01C Earth 
gravity model 5 |25j . the systematic error due to the remaining uncancelled 
even zonal harmonics J4, Jq, Jg... amounts to 5%. This result has been 
obtained with a 1-a RSS (Root Sum Square) calculation based on the vari- 
ance matrix of EIGEN-CG01C because, contrary to EGM96 and the other 
previous Earth gravity models, in this case the recovered multipolar coeffi- 
cients are well disentangled. A pessimistic 1-a upper bound of 6% can be 
obtained by summing the absolute values of the individual errors. 

Note also that the combination © preserves one of the most important 
features of the combination (J2J) of orbital residuals: indeed, it allows to cancel 
out the very insidious 18.6-year tidal perturbation which is a t = 2, m = 

4 The possibility of using only the nodes of the LAGEOS satellites in view of the more 
accurate GRACE Earth gravity solutions was presented for the first time in |16|. although 
without quantitative details. 

5 Such model represents a long-term averaged solution which combines data from 
CHAMP (860 days), GRACE (109 days) and surface gravimetry/altimetry; moreover, 
the released sigmas of the spherical harmonic coefficients of the geopotential are not the 
mere formal statistical errors, but are calibrated, although preliminarily. Then, guesses 
of the impact of the systematic error due to the geopotential on the measurement of the 
Lense-Thirring effect based on this solution should be rather realistic. However, caution 
is advised in considering the so obtained evaluations because of the uncertainties of the 
calibration process which affect especially the even zonal coefficients |17|. 
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constituent with a period of 18.6 years due to the Moon's node and nominal 
amplitudes of the order of 10 3 mas on the nodes of LAGEOS and LAGEOS II 
|26| . On the other hand, the impact of the non-gravitational perturbations 
on the combination (JSJ) over a time span of, say, 7 years could be quantified 
in 0.1 mas yr _1 , yielding a 0.3% percent error. The results of Table 2 and 
Table 3 in have been used. It is also important to notice that, thanks to 
the fact that the periods of many gravitational and non-gravitational time- 
dependent perturbations acting on the nodes of the LAGEOS satellites are 
rather short, a reanalysis of the LAGEOS and LAGEOS II data over not 
too many years could be performed. As already pointed out, this is not 
so for the combination (J2J) because some of the gravitational j2@] and non- 
gravitational |21| perturbations affecting the perigee of LAGEOS II have 
periods of many years. 

A possible weak point of (j5J) could be represented by the fact that it is 
affected by the secular variations of the uncancelled even zonal harmonics 
J4, Jq, .... For the topic of Jj>, which has recently received great attention by 
the geodesists' community in view of unexpected variations of 6 J2, see [2"7] , 
The mismodelled shift of ® due to the secular variations of the uncancelled 
even zonal harmonics can be written as 

£ (^ AGEOS + ft ^ AGEOS n ) ^T 2 bs , (6) 
1=2 

where the coefficients Qj are <90 c i ass /<9J£ and have been explicitly calculated 
up to degree t = 20 in ^5]. It must be divided by the gravitomagnetic shift 
of © over the same observational time span 

/qLAGEOS , _ qLAGEOS iA rp /l 7 o mQO ,,,.-1 t (7\ 

I ' 'lt +<Zi i 'LT I Jobs — 47.9 mas yi i ODS - (7) 

By assuming 5J4 = 0.6 x 10~ n yr -1 and SJq = 0.5 x 10~ n yr" 1 (Cox et al 
2002), it turns out that the 1-a percent error on the combination which 
grows linearly with T OD s, would amount to 1% over one year. However, it 
must be pointed out that it is very difficult to have reliable evaluations of the 
secular variations of the higher degree even zonal harmonics of geopotential 
also because very long time series from the various existing laser-ranged 
targets are required. 

Recently, ((SJ has been adopted for a 11-years analysis j^H] based on the 
2nd generation GRACE-only EIGEN-GRACE02S Earth gravity model [23] . 

6 Fortunately, any issues concerning J2 do not affect the combination JjJJ. 
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although the claimed total accuracy of the performed test might be rather 
optimistic |3TJ| . 

Since GRACE is expected to improve especially the mid to high degree 
part of the spectrum of the even zonals, in [2U] a J2 — ... Js— free combination 
involving the nodes of the spherical geodetic satellites LAGEOS, LAGEOS 
II, Ajisai, Starlette and Stella 7 is proposed. This combination would be 
affected just by the higher degree Je and would be insensitive to the secular 
variations of the lower degree even zonal harmonics. However, due to the 
low altitude of Starlette and Stella, the 1-a error due to the static part of 
the geopotential would still amount to 10% (RSS calculation) with a 1-a 
upper bound of 31%, according to the present-day EIGEN-CG01C model. 
It may be interesting to note that the combination (J2J) would be affected 
less than 1% by the mismodelling in the static part of geopotential. 

4 The use of the laser- ranged satellite Ajisai and 
of the radar altimeter satellite Jason-1 

Up to now, the attention has been focused exclusively on the spherical laser- 
ranged satellites because of the high accuracy with which it is possible to 
determine their orbits. It is so also thanks to their high altitude, small 
area-to-mass S/M ratio and spherical shape, which reduces the impact of 
the non-gravitational perturbations. A recent, important achievement in or- 
bit determination regards Jason-1. It is a radar altimeter satellite, launched 
on December 7, 2001, as a follow-on to the very successful TOPEX/Poseidon 
mission which was launched in 1992. The satellite carries state-of-the-art 
hardware for the three most accurate tracking systems available: Satellite 
Laser Ranging (SLR), Doppler Orbitography and Radiopositioning Inte- 
grated by Satellite (DORIS) and Global Positioning System (GPS). In ad- 
dition, its radar altimeter measurements over the oceans can be used for 
either independent orbit validation, or as additional tracking data. In it 
has been shown that it is possible to reach the 1-cm level in determining its 
orbit (in the radial direction) by using this dense coverage of precise tracking 
data. 

This is a very important result, because the orbital parameters of Jason- 
l 8 make it suitable to be used, in principle, for the purpose of an accurate and 

7 The current RMS of fit of the computed orbits to the laser data of these satellites is 
of the order of 1 cm or even better, especially for the LAGEOS satellites. 

8 One could argue that TOPEX/Poseidon, with the same orbital parameters as Jason- 
1, a smaller area-to-mass ratio and a 12-year (and counting) lifetime, would be the more 
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reliable measurement of the Lense-Thirring effect, provided that accurate 
force models are available for a dynamic orbit determination. Indeed, the 
Jason-1 orbit parameters are rather similar to those of the geodetic Ajisai 
satellite. In Tabled the orbital parameters of the LAGEOS satellites and of 
Ajisai and Jason-1 are reported. However, it must be preliminarily pointed 
out that special care should be taken in handling the orbital manoeuvre 
burns of Jason-1, which are designed to counteract the drifting of the orbit 
and keep it on its repeating ground-track, as well as the infrequent safe- 
mode periods, which also complicate the dynamic modelling. These events 
do not prevent us from measuring any long-term drifts in the node however, 
as long as the orbital arcs are designed to start after and stop just before 
these instances. Because the manoeuvres are mostly performed in pairs, 
approximately one hour apart, and the two safe-mode periods up to now 
have lasted several days, they would however introduce uncertainties when 
fitting for long-period signals. 

Table 1: Orbital parameters, predicted Lense-Thirring nodal rates Olt an d 
approximate area-to-mass ratios S/M. a, e, i and n are the semimajor axis, 
the eccentricity, the inclination to the Earth's equator and the Keplerian 
mean motion n = \J GM/a 3 , respectively. 



Satellite 


a (km) 


e 


i (dcg) 


n (1(T 4 s- 1 ) LT (mas yr- 1 ) 


S/M (m 2 kg- 1 ) 


LAGEOS 


12270 


0.0045 


109.84 


4.643 30.7 


6.9 xHT 4 


LAGEOS II 


12163 


0.0135 


52.64 


4.710 31.4 


7.0 xl0~ 4 


Ajisai 


7870 


0.001 


50.0 


9.042 116.2 


5.3 xlO~ 3 


Jason-1 


7713 


0.0001 


66.04 


9.320 123.4 


2.7 xl0~ 2 



4.1 A possible combination of nodes and the gravitational 
errors 

The following combination 

^LAGEOS + ^LAGEOS II + ^Ajisai + ^Jason-1 = ^49^ (g) 
with 

fci = 0.347, k 2 = -0.005, k 3 = 0.068, (9) 

suitable satellite. However, TOPEX/Poseidon also has a more irregular attitude and solar 
array behavior, a more complex shape and less precise tracking instrumentation. Its RMS 
radial orbit accuracy is estimated at 2.8 cm |33) . 
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could be used in order to cancel the effect of the first three even zonal 
harmonics J2, J4 and Jq along with their temporal seasonal, stochastic and 
secular variations. Recently, an analogous approach has been followed |32| 
by including also some GPS satellites. According to EIGEN-CG01C, the 
systematic error induced by the remaining uncancelled even zonal harmonics 
of geopotential amounts to 0.7% only (RSS calculation at 1-a level), with 
an upper bound of 1.6% (sum of the absolute values of the individual errors 
at 1-a level). It must be noted that, contrary to the combination which 
can be obtained with the nodes of Stella and Starlette instead of that of 
Jason-1, a calculation up to degree £ = 20 is quite adequate to assess the 
error due to geopotential. It is so thanks to the higher altitude of Jason-1 
with respect to Stella (asti = 7193 km) and Starlette (astr = 7331 km). 
Contrary to the other examined combinations, in particular the LAGEOS- 
LAGEOS II two-nodes J2-free formula of [20], the secular variations of the 
even zonal harmonics Je do not affect (JHJ). The combination Q would not 
be affected by harmonic perturbations of tidal origin with large amplitudes 
and particularly long periods. Indeed, the tidal perturbation induced by the 
tesseral (m = 1) K\ tide, which is one of the most powerful tidal constituents 
in affecting the orbits of Earth satellites has its period equal to that 
of the satellite's node; it amounts to -0.47, -0.32, -1.55 and 2.84 years for 
Jason-1, Ajisai, LAGEOS II and LAGEOS, respectively. This means that 
a relatively short observational time span could be adopted in order to have 
it longer than the periods of the most effective perturbations which could 
then be fitted and removed from the time series |34| I35j . 

4.2 The impact of the observational errors of Ajisai and 
Jason-1 

In regard to the impact of the measurement errors, let us assume, in a very 
pessimistic way, that the RMS of the recovered orbits of Ajisai and Jason-1, 
obtained in a truly dynamical way without resorting to too many empirical 
accelerations which could absorb also the Lense-Thirring effect, amount to 
1 m over, say, 1 year. Then, the error in the nodal rates can be quantified 
as 26.2 mas and 26.7 mas for Ajisai and Jason-1, respectively. Thanks to 
the coefficients ©, their impact on the combination (JHJ) would amount to 
1.6 mas, i.e. 3.4% of the Lense-Thirring effect over 1 year. 
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4.3 The impact of the non-gravitational perturbations of Aji- 
sai 

Ajisai is a spherical geodetic satellite launched in 1986. It is a hollow sphere 
covered with 1436 corner cube reflectors (CCR's) for SLR and 318 mirrors 
to reflect sunlight. Its diameter is 2.15 m, contrary to LAGEOS which has 
a diameter of 60 cm only. Its mass is 685 kg, while LAGEOS mass is 406 
kg. Then, the Ajisai's area-to-mass ratio S/M, which the non-conservative 
effects are proportional to, is larger than that of the LAGEOS satellites by 
almost one order of magnitude resulting in a higher sensitivity to surface 
forces. However, we will show that their impact on the proposed combination 
© should be less than 1%. 

4.3.1 The atmospheric drag 

The most important non-conservative force in affecting the orbits of the low 
satellites is the atmospheric drag. Its acceleration can be written as 

«o = -lcn^)pVV, (10) 

where Co is a dimensionless drag coefficient close to 2, p is the atmospheric 
density and V is the velocity of the satellite relative to the atmosphere 
(called ambient velocity). Let us write V=v— u>x r where v is the satellite's 
velocity in an inertial frame. If the atmosphere corotates with the Earth 
er is the Earth's angular velocity vector = u<^k, where k is a unit 
vector. However, it must be considered that there is a 20% uncertainty in 
the corotation of the Earth's atmosphere at the Ajisai's altitude 9 . We will 
then assume a = w®(l+£)fc, with £ = 0.2 in order to account for this effect. 

Regarding the impact of a perturbing acceleration on the orbital motion, 
the Gaussian perturbative equation for the nodal rate is 

-jr = 7= s . : a n 1 - J sum, (11) 

«r nay 1 — e z sin % 

where n = a/ GM/a 3 is the Keplerian mean motion, An is the out-of-plane 
component of the perturbing acceleration and u = co + / is the satellite's 
argument of latitude. 

9 It is believed that the atmosphere rotates slightly faster than the Earth at some 
altitudes with a 10-20% uncertainty (Ries 2004, private communication) . 
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The out-of-plane acceleration induced by the atmospheric drag can be 
written as [3*5] 

^(atm) _ _ J^ D(J p Vr s [ n j^ cosu ^ (12) 



with 



2ah cos i (or cos 5 x 



and 

kR = 1 _ ~- p- + | 11J2 ^ | _ (M) 
f z \ w / 

The quantities /i and 5 are the orbital angular momentum per unit mass 
and the satellite's declination, respectively. By inserting (fT2|) in (fTT|) and 
evaluating it on an unperturbed Keplerian ellipse it can be obtained 

jq (atm) -i 

^ K--K DP {l-e*)oa. (15) 

It must be pointed out that the density of the atmosphere p has many 
irregular and complex variations both in position and time. It is largely 
affected by solar activity and by the heating or cooling of the atmosphere. 
Moreover, it is not actually spherically symmetric but tends to be oblate. 
A very cumbersome analytic expansion of p based on the TD88 model can 
be found in [HE]- I n order to get an order of magnitude estimate we will 
consider a typical value p = 1 X 10~ 18 g cnr 3 at Ajisai altitude [STj. By 
assuming Cd = 2.5 ()15j) yields a nominal amplitude of 25 mas yr _1 for 
the atmosphere corotation case and 5 mas yr _1 for the 20% departure from 
exact corotation. The impact of such an effect on Q would be 3 x 10~ 3 . 



4.3.2 The thermal and radiative forces 

The action of the thermal forces due to the interaction of solar and terrestrial 
electromagnetic radiation with the complex physical structure of Ajisai has 
been investigated in 3?.. The temperature asymmetry on Ajisai caused 
by the infrared radiation of the Earth produces a force along the satellite 
spin axis direction called Yarkovsky-Rubincam effect. This thermal thrust 
produces secular perturbations in the orbital elements, but no long-periodic 
perturbations exist if the spin axis of Ajisai is aligned with the Earth's 
rotation axis. In fact, the spin axis was set parallel to the Earth rotation axis 
at orbit insertion. The analogous solar heating (Yarkovsky-Schach effect) 
is smaller than the terrestrial heating. A nominal secular nodal rate of 15 
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mas yr 1 due to the Earth heating has been found. It would affect ([8j) at a 
1.5 x 10" 3 level. 

The effect of the direct solar radiation pressure on Ajisai has been studied 
in |38| . For an axially symmetric, but not spherically symmetric, satellite 
like Ajisai there is a component of the radiation pressure acceleration di- 
rected along the sun-satellite direction a^°^ and another smaller component 
perpendicular to the sun-satellite direction. By assuming for the isotropic 
reflectivity coefficient its maximum value Cr = 1.035 the radiation pressure 
acceleration Oq°^ experienced by Ajisai amounts to 2.5 x 10~ 8 m s~ 2 . Its 

nominal impact on the node, proportional to ea^°^ /4na, can be quantified 
as 5.7 mas yr -1 ; it yields a 5.7 x 10 -4 relative error on (|8}l. The anisotropic 
component of the acceleration would amount, at most, to 2% of the isotropic 
one, so that its impact on (JHJ) would be totally negligible. 

4.4 The impact of the non-gravitational perturbations on 
Jason-1 

The complex shape, varying attitude modes and the relatively high area-to- 
mass ratio of the satellite (see Table ^) , suggest a more complex modelling 
and higher sensitivity to the non-gravitational accelerations than in the case 
of the spherical geodetic satellites. On the other hand, important limiting 
factors in the non-gravitational force modelling of the spherical satellites 
are attitude and temperature knowledge |21| . These parameters are actually 
very well-defined and accurately measured [39 j on satellites such as Jason-1. 
A lot of effort has already been put into the modelling of non-gravitational 
accelerations for TOPEX/Poseidon [301 EH EH , so that similar models [23| 
have been routinely implemented for Jason-1. These so-called box-wing 
models, in which the satellite is represented by eight flat panels, were de- 
veloped for adequate accuracy while requiring minimal computational re- 
sources. A recent development is the work on much more detailed models 
of satellite geometry, surface properties, eclipse conditions and the Earth's 
radiation pressure environment for use in orbit processing software |44l ?] . It 
should be noted that such detailed models were not yet adopted in the orbit 
analyses of |31j . In fact, their results were based on the estimation of many 
empirical 1-cycle-per-revolution (cpr) along-track and cross-track acceler- 
ation parameters, which absorb all the mismodelled/unmodelled physical 
effects, of gravitational and non-gravitational origin, which induce secular 
and long-period changes in the orbital elements. Due to the power of this 
reduced-dynamic technique, based on the dense tracking data, further im- 
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Figure 2: Time series of the modelled Jason-1 out-of-plane solar radiation 
pressure acceleration (grey) and the estimated 1-cpr out-of-plane accelera- 
tion (black, right-hand scale). 



provements in the force models become largely irrelevant for the accuracy 
of the final orbit. Such improved models remain, however, of the highest 
importance for the determination of <5f2 Jason_1 . 

From it can be noted that, since we are interested in the effects av- 
eraged over one orbital revolution, the impact of every acceleration constant 
over such a timescale would be averaged out. The major problems come from 
1-cpr out-of-plane accelerations of the form = Sn sin u + CV cos u, with 
Sn and CV constant over one orbital revolution. 

We have analyzed both the output of the non-gravitational force mod- 
els described in and the resulting residual empirical 1-cpr accelerations 
estimated from DORIS and SLR tracking over 24-hour intervals. Solar 
radiation pressure, plotted in figure |2J is by far the largest out-of-plane 
non-gravitational acceleration, with a maximum amplitude of 147 nm s~ 2 . 
It is followed by Earth radiation pressure at approximately 7 nm s~ 2 . The 
contributions of aerodynamic drag and the thermal imbalance force on the 
cross-track component are both estimated to have a maximum of approxi- 
mately 0.5 nm s~ 2 . As can be seen in figure|2j the cross-track solar radiation 
pressure acceleration shows a sinusoidal long-term behavior, crossing zero 
when the Sun-satellite vector is in the orbital plane, roughly every 60 days. 
It is modulated by the long-term seasonal variations in Sun-Earth geom- 
etry, as well as by eclipses and the changing satellite frontal area, both of 
which are contributing 1-cpr variations. In fact, the shading of certain areas 
in figure [21 is due to the effect of the eclipses, which, at once-per-orbit, occur 
much more frequently than can be resolved in this figure. 

As mentioned before, the empirical 1-cpr accelerations absorb the errors 
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of almost all the unmodelled or mismodelled forces. Now note the systematic 
way in figure [21 in which the empirical 1-cpr cross-track acceleration drops 
to values of below 1 nm s~ 2 near the end of each eclipse-free period, and has 
its maximum level of 5-6 nm s~ 2 only during periods containing eclipses. 
The fact that the amplitude, but also the phasing (not shown in the figure) 
of the 1-cpr accelerations show a correlation with the orientation of the 
orbital plane with respect to the Sun, indicates that it is for a large part 
absorbing mismodelled radiation pressure accelerations. 

By averaging (jllj) over one orbital revolution and from the orbital pa- 
rameters of Table ^ it turns out that a 1-cpr cross-track acceleration would 
induce a secular rate on the node of Jason proportional to 7.6 x 1CP 5 s 
m _1 x Sn m s~ 2 . This figure must be multiplied by the combination coeffi- 
cient k^. By using the average value of the empirical 1-cpr acceleration from 
the above analysis Sn ~ 2.3 nm s~ 2 as an estimate for the mismodelled non- 
gravitational forces, it can be argued that the impact on fc3<5f2 Jason ~ 1 would 
amount to 77.4 mas yr _1 . However, it must be pointed out that our assumed 
value of SV can be improved by adopting the aforementioned more detailed 
force models or by tuning the radiation pressure models using tracking data. 
In addition, it must be pointed out that Sn experiences long-term varia- 
tions mainly induced by the orientation of the orbital plane with respect to 
the Sun, and the related variations in satellite attitude. For Jason-1 such a 
periodicity amounts to approximately 120 days (the (3 cycle). Let us, now, 
evaluate what would be the impact of such a long-periodic perturbation on 
our proposed measurement of the Lense-Thirring effect. Let us write, e.g., 
a sinusoidal law for the long-periodic component of the weighted nodal rate 
of Jason-1 




K3— — = (77.4 mas yr J x cos 
at 

then, if we integrate (j!6|) over a certain time span T \, s we get 



/, : .,A<> ( y~ ) (77.4 mas yr -1 ) x sin 



2tt 



(16) 



(17) 



Then, the amplitude of the shift due to the weighted node of Jason-1, by 
assuming Pn< =120 days, would amount to 

k 3 Att < 4 mas. (18) 

The maximum value would be obtained for 

■^ = 7, 3 = 1,3, 5, ...5*30,90,150,... days. (19) 
4 
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So, the impact on the proposed measurement of the Lense-Thirring effect 
would amount to 



^LT S RP 



< 



(4 mas) 



(20) 



(49.5 mas yr x ) x (T obs yr) 



for, say, T Q b s = 2 years (f2*0|) yields an upper bound of 4%. Moreover, it must 
also be noted that it would be possible to fit and remove such long-periodic 
signals from the time series provided that an observational time span longer 
than the period of the perturbation is adopted, as done in [TT] . 

5 Conclusions 

According to the EIGEN-CG01C Earth gravity model, which combines data 
from CHAMP, GRACE and surface gravimetry/altimetry, the systematic 
error in the measurement of the Lense-Thirring secular effect due to the 
mismodelling in the static part of the geopotential is close to a few per- 
cent level, at 1-a level, if the combinations (JHJ) is to be used. From this 
point of view, ©, which is affected by Js, Jio--- at 1.6% level (1-<t), is, at 
present, better than © which, instead, is affected by J4, Jq, Js, ... at a 6% 
level (1-c). Another advantage of © with respect to © is that the com- 
bination involving Jason- 1 and Ajisai is insensitive to the secular variations 
of the even zonal harmonics, contrary to the two-nodes LAGEOS-LAGEOS 
II combination which is, instead, affected by the mismodelling in J4 and Jq. 
This could represent a drawback because their effect grows linearly in time 
and, at present, the knowledge of Ji is rather poor. Their estimated impact 
on (J5J) amounts to 1%, at 1— a level, over one year; time spans of more than 
just one year are required in this kind of analyses. 

The price to be paid for the benefits in the reduction of the systematic 
error of gravitational origin is mainly represented by the non-gravitational 
perturbations introduced along with the new orbital elements. However, 
while the bias due to the geopotential is unavoidable because it is constant 
or grows linearly in time, the signature of the non-conservative accelerations 
is long-periodic, so that they could be fitted and removed from the time se- 
ries. This perspective should strongly encourage to explore the possibility 
of a reanalysis of the Jason-1 orbital data in order to make them suitable 
to be combined with those from the spherical geodetic satellites. Particular 
attention should be paid to an as accurate as possible truly dynamical mod- 
elling of the non-gravitational accelerations acting on the node of Jason-1. 
This would be helpful also in other scientific branches related to the Jason-1 
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activity. It must be noted the non-gravitational perturbations, which are 
almost negligible on the nodes of the LAGEOS satellites, could represent a 
major obstacle in using Jason- 1. A reasonable estimate of the impact of the 
non-gravitational perturbations modelling error on Jason-1 on the combi- 
nation (jHJ) yields an error of the order of 4% over an observational time span 
of 2 years. However, since the main periodicity of the direct solar radiation 
pressure perturbation seems to amount to approximately 120 days only, over 
a time span of a few years it should be possible to average it out or, at least, 
to fit and remove such biasing signature from the time series. Of course, 
if the combination (|SJ) will be used, a careful choice of the observational 
time span would be required in order to reduce the uncertainties related to 
the orbital maneuvers, which, however, are mainly-although not entirely-in 
plane, and safe-mode events of Jason-1. 

Finally, the obtainable accuracy with the node-node-perigee combina- 
tion (J2j), whose error due to geopotential will remain smaller than that of 
is strongly related to improvements in the evaluation of the non- 
gravitational part of the error budget and to the use of time spans many 
years long. However, it neither seems plausible that the error due to the 
non-conservative forces acting on the perigee of LAGEOS II will fall to the 
1% level nor that a reliable and undisputable assessment of it will be eas- 
ily obtained. Alternative combinations including the orbital data from the 
other existing SLR geodetic spherical satellites are, at present, not compet- 
itive with the combinations (J5J) and (JHJ); indeed, the systematic error due to 
the static part of geopotential amounts to 31% at 1— a level. 
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